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Abstract: Twisted Eguchi-Kawai reduced chiral models are shown to be formally 
equivalent to a U{1) non-commutative parent theory. The non-commutative theory 
describes the vacuum dynamics of the non-commutative charged tachyonic field of 
a brane system. To make contact with the continuum non-commutative theory, a 
double scaling large N limit for the reduced model is required. We show a possible 
limiting procedure, which we propose to investigate numerically. Our numerical 
results show substantial consistency with the outlined procedure. 
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1. Introduction 

Twisted Eguchi-Kawai (TEK) reduced models ^ provide a non-perturbative defi- 
nition of certain non-commutative field theories p, ^ It has been shown that 
the U (N) gauge theory with gauge fields obeying twisted boundary conditions over 
the non- commutative torus Tq is equivalent to a U{N) gauge theory, with N suit- 
ably chosen, over the non-commutative torus Tq, with gauge fields obeying periodic 
boundary conditions [Q. This is a consequence of a more general fact, known as 
Morita equivalence. 

We would like to propose an application of the formalism of ref. Q to principal 
chiral models, possibly providing numerical evidence, in order to show that two- 
dimensional TEK reduced chiral models can be considered as a non-perturbative 
definition of a non-commutative field theory. 

Explicitly, we are going to show the equivalence of TEK principal chiral model 
with symmetry group U{N) to a non-commutative U{1) lattice theory compactified 
on a torus with periodic boundary conditions. We will describe the correspond- 
ing non-commutative theory, i.e. its action and symmetries, coupling constant and 
dimensionful non-commutativity parameter. We will eventually try to define a pro- 
cedure which can lead to a sensible continuum limit, and therefore define the set of 
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values of the coupling constants and of the symmetry parameter N of the original 
TEK theory relevant for a numerical check of the consistency of the whole approach. 

The reduced model, in this context, should be considered in a different limit from 
the original planar one. In the planar limit one has to send N ^ oo while keeping 
the lattice coupling constant f3 fixed. The continuum limit is then reached as one 
sends (3 ^ oo, because the RG relation a oc e~'^'^ exists between the lattice spacing 
a and the lattice coupling constant (3. The limit is called planar since the large N 
dominant Feynman graphs are, in this limit, planar, i.e. they can be drawn on a 
Riemann surface of genus g = 0. 

For the TEK model to reproduce the corresponding non- commutative field the- 
ory we must consider a different limiting procedure to approach the continuum limit. 
This procedure is called double scaling limit, since (3 and N must be sent to infinity 
in a correlated manner. We know that in the non reduced original field theory this 
corresponds to taking into consideration also non-planar graphs, analogously to what 
happens in matrix models of 2D gravity: contributions from higher genus topologies 
imply a higher symmetry in the problem, which enables to make contact with the the- 
ory one wants to reproduce. Quite the same thing happens in non-commutative field 
theory, since the interaction vertex is invariant only up to cyclic permutations of the 
momenta, and therefore one needs to keep track of the cyclic order in which lines em- 
anate from vertices in a given Feynman diagram. Non-commutative Feynman graphs 
are thus ribbon graphs drawn on Riemann surfaces of particular genus, in complete 
analogy with what happens in the ordinary large N limit of field theories. Non- 
planar contributions thus naturally arise in the present context of non-commutative 
field theories. 

Principal chiral models in D = 2 are since long known to be in a sense^ 
a simpler counterpart of lattice Yang-Mills theory in D = 4. Also as far as the 
reduced TEK models of the two theories are concerned, the case of chiral models 
is more tractable, and even numerical simulations are easier and more conclusive. 
For this reason we consider principal chiral model as an important test for the non- 
commutative interpretation of TEK reduced models. 

Furthermore, the non-commutative theory arising from the TEK reduced prin- 
cipal chiral model turns out to possess exactly the action needed to describe the 
vacuum dynamics of the non-commutative charged tachion of a certain brane system 

i- 

We begin with a general description of non- commutative field theories, in order 
to fix the notation, and describe the particular non-commutative field theory we 
are interested in, i.e. 2D principal chiral models. We then pass to a short analysis 
of the reduced models and of the consequences of our numerical results. In the 
final section we demonstrate the equivalence between TEK reduced chiral models 

^This correspondence is exact in D = 1 and D = 2 respectively. 
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and non-commutative U{1) principal chiral models on the lattice, indicating the 
procedure needed to approach a sensible continuum limit. This will eventually lead 
to a concrete proposal of a numerical study of the theory reread in this new light. 
The physical interpretation of the presented non-commutative theory is then briefly 
sketched. We include the results of our Monte Carlo simulations, which directly 
enforce the correctness of the whole procedure. 

2. Non-commutative Field Theories 

We will briefly recall the so called Weyl quantization procedure for field theories 
on non-commutative spaces. Let us consider for simplicity a scalar field theory on 
Euclidean M^, defined by some action S[(j)], and whose partition function is as usual 

-5[0] 



V(l) e-^'^^K (2.1) 

In order to pass from ordinary to non- commutative space-time, we replace the local 
coordinates by hermitian operators a;^ which have the following commutations 
relations: 

[x^,x„] = iOf,^, (2.2) 

where 6^i, = —9^^ is a real valued anti-symmetric matrix with the dimensions of 
length squared. Implementing the substitution x^ x^ we obtain the Weyl operator 
(f) corresponding to the field 0. 

One can define an operator which transforms a field into its Weyl operator: 

(f) = J d^x(f){x)A{x), 

(2.3) 



It is possible to introduce an anti- hermitian derivation 5^ through the commutations 
relations 

[d^,x,]=6^,, [d^,d,]=0. (2.4) 
The operator obviously satisfies the following relations: 

[d,J] = Jd^xd,cf>{x)A{x), 

(2.5) 

[a^,A(x)] = d^A{x). 
This implies that, given the generator of the translations e''''^'', G M, which satisfies 

e"'^^' A(x)e-''^^'' = A{x + v), (2.6) 
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the operation TrA(a;) is independent of x for any trace on the algebra of operators. 
Therefore the integration of the fields on the space-time is represented by 

Tr0 = j d^x <f){x), TrA{x) = 1. (2.7) 
This enables to define the inverse of the correspondence between fields and operators: 

(f){x) = Tr (4)A{x)^ . (2.8) 

In conclusion, in order to pass from ordinary to non-commutative field theory, we 
need to implement the following substitutions: 

0^0, d,<P-^[d„4>]. (2.9) 

The substantial difference between the two theories comes from the definition of 
products of fields. In the non-commutative case one has, for 4>3 = 4>i4^2 



(l)3{x) = IV (0i02A(x)) = 

= t/d^^0i(y)02(^)e-2^(^"')'"'(^-^)''(^-^)-' = 

(2.10) 

= (f)i{x) exp (^^%9^^dl^ (f)2{x) = (pi{x) -k (t)2{x) = 

= ct>,{x)Ux) + Er=i (I)" if^^"^ • • • ^'"'"^n • • • ■ ■ ■ a,>2(x), 

which defines the star or Moyal product of the fields. 

2.1 Non-commutative U{N) Principal Chiral Models 

The theory we would like to study is defined in a D dimensional Euclidean space by 
the following action, multiplied by a suitable coupling constant 

S = (3N J d^x tr(^r) {d^U{x)d''U\x)) , (2.11) 

where the sum over = 1,2, ...,D, and a flat Euclidean metric are intended^. The 
unitary matrices Uij{x), with i,j — 1, 2, N satisfy the condition 

{U\x)U{x))^^ = 5., Vx. (2.12) 
The partition function is defined over the usual Haar measure 

dU e-^[^l. (2.13) 



^In what follows we will restrict to the case D = 2. 
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The action is naturally invariant under the transformations (where Vl, Vr G U{N)) 



U VlU, U ^ UVr. 
In the present formalism, one needs to define the Weyl operators as 



(2.14) 



11 = j d^x A{x) ■ U{x), (2.15) 
where the operator A[x) is the one defined in ( p.3|) . The action is rewritten as 



S = (3NTt[ tr(;v) ( 



(2.16) 



where Tr is the trace operator over coordinates, while tr(N) is the (finite-dimensional) 
trace in the fundamental representation of the U (N) group. 

The corresponding non-commutative fields U{x), defined through the inverse trans- 
formation ( p.8|) , satisfy the star-unitarity condition^ 



U{x) ^ U{xy = U{x)^ * U{x) = In. 
The action is given by the obvious translation of eq. ( p.l6| ) 



(2.17) 



S = 13N j d^x tr(jv) {d^U{x)i<d^'U{xy) . (2.18) 
The invariance of the action naturally reads, owing to cyclicity of the trace tr(7v), 

U{x) ^ ql U{x) U{x) ^ U{x) qr, (2.19) 

where the N x N matrices QljQr are ordinary unitary matrices, therefore satisfying 
99^ = ^N- 



3. The reduced model 

The principal chiral model, in Z) = 2, is defined on the lattice via the usual substi- 
tution of the derivative with a finite difference 



d,U{x) 



(3.1) 



and the resulting action reads 



S = -PNJ2J2 + Ux^,Ul\ . (3.2) 

X fj.=l,2 

The naive Eguchi-Kawai reduction prescription — > t/ is clearly not applicable in 
this context. Instead, one can resort to the TEK prescription, which is defined as 

^since det{f -k g) ^ dct(/) -kdet{g) we cannot take SU{N) as the symmetry group of the theory. 
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^ T^TTUiTlriTlr, (3.3) 
where the twist matrices obey the Weyl-'t Hooft algebra 

r.r^, (3.4) 

where is the parameter of the symmetry group U{N) or SU{N) (and thus of the 
matrix U) and n^^/ is an integer valued antisymmetric tensor, whose generic form in 
D = 2 is of course 



exp 



27ii 



= ( °^ ^ ) , MeZ. (3.5) 



For a given and M the solution to ( p.4|) is provided, up to global SU{N) trans 



formations, by the N x N shift and clock matrices^ 

MM) _ r- 

The two matrices will be given respectively by S and C. 
Applying the reduction prescription to the action (|3.2D gives 



(3.6) 



Stek = -PN ^ Tr [UT.U^TI + /i.e.] . (3.7) 
We notice that the model shows two symmetries, namely: 



2. U ^ z-U, z el. 



The first symmetry is reminiscent of the space-time translational symmetry of the 
original model (indeed from ( |3.3|) it is clear that the role of the F^ is that of generators 
of translations in the dual lattice of the reduced theory) while the second represents 
the residual global symmetry SU (N) x SU{N) of the parent theory reduced to the 
center Z^r of the algebra of SU{N). In the case of a symmetry U{N) x U{N) the 
second symmetry is of course a U{1) symmetry. 

The introduction of the TEK reduced model was originally motivated by its 
supposed equivalence with the parent theory (defined by the action of eq. (|3.2|) ) in 
the large A^ limit. This equivalence should follow basically from two facts: 



^The two matrices are a natural extension of the known 't Hooft twist matrices 
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- The Schwinger-Dyson (SD) equations of the reduced theory and that of the 
parent theory are exactly the same in the large limit, up to terms which are 
not invariant under the Z^y symmetry 

- If in all regimes the symmetry is not spontaneously broken, the two models pos- 
sess exactly the same SD set of equations, and given the same initial conditions 
should coincide 

It has been found [^j that the invoked equivalence holds for the strong coupling region, 
while in the weak coupling the reduced theory is manifestly different from its parent 
version. Thus, although the second symmetry does not seem to be spontaneously 
broken, the SD based argument for the equivalence is not sufficient, at least in a 
certain regime, to completely identify the theory. 

It therefore makes sense to ask whether the reduced principal chiral model can 
provide a non perturbative definition of another theory. Our claim is that this 
theory can be interpreted as the non- commutative U{1) version of ( |3.2| ), and that, 
in a suitable limit to be defined, it reproduces its continuum limit, i.e. the one 



specified in the Weyl operator notation by eq. ( p.l6|) and in the non-commutative 
field notation by eq. (|2.18|) . 



4. Reduced chiral models and non-commutative chiral fields 

We will show in what follows how to map the TEK model into the non-commutative 
version of the theory defined by the action (|3.2|) . Given an integer valued vector 
k = [ki, k2), we introduce the N x N matrices 

_ ■pA:ipA:2gVr«(ni2)A:2fci/A'' _ -pA;i-pfc2g7riMfe2'si/-'V 

The phase factor is given to symmetrically order the product of twist eaters. Inci- 
dentally, the Jfc's have the same algebraic properties as the plane Weyl basis e*'^*'^' 
for the continuum non- commutative field theory on the torus [Q. 
The relevant properties of these matrices are that there are only N"^ such matrices, 
owing to the periodicity properties 

JN-k = J-k = JI, (4.2) 

and that they obey the orthonormality and completeness relations 

Jft'!C(N) i^JkJ\) = 5/c,g(mod N), 



jf^k€Z%i'^k)tiu{Jk)xp — S/ipSuX- 



(4.3) 



They therefore form a basis for the linear space gl{N, C) of N x N complex matrices, 
and in particular one can expand a matrix U as 
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^ = iV^ E f^(^)'^'^' f^(^) = ^^^W i^-Jl) ■ (4.4) 

We can interpret the momentum coefficients as the dynamical degrees of freedom in 
the TEK model 

In analogy with the continuum counterpart ( p.3|) we can define the operator (this 
time on a discrete torus) 

where L = aN is the dimensionful extension of the lattice with A^^ sites x*. Because 
of the relations ( [4.2p the A{x) matrices are Hermitian and periodic in with period 
L, and thus the lattice is a discrete torus. 

In analogy with the continuum formalism depicted in section ^, we can define 
an invertible map between N x N matrices and lattice fields. Namely, we have the 
following properties 

tr(^) (JfcA(x)) = ^e^-^^'-Vi, 

A(x)^,A(x)ap = 6^p6,x, (4.6) 

^tr(Ar) (A(x)A(?/)) = A^^(5^,j,(mod L)- 

which yield a natural definition for the lattice field U{x) from the Fourier modes of 
its matrix partner U: 

^(^) = ^ E me'""""^'/' = ^trw {UA{x)) . (4.7) 

Since 

X 

the unitarity condition on the matrix U is translated on the field U{x) in terms of 
U{1) star unitarity: 

U{x)^U*{x) =U*{x)^l{{x) = 1, (4.9) 

where the lattice star product is defined by the natural discrete analog of eq. ( |2.1CI| ), 
namely 

A^B = j^tnM) {ABA{x)) 

(4.10) 
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with dimensionful non-commutativity parameter 



0^,, = n^,. (4.11) 

vr 



The star product ( [4. 101) reproduces the continuum version of eq. (|2.10| ) in the hmit 



a — i> 0, and it reproduces the same algebraic properties with space-time integrals 
replaced by lattice sums. 

In order to reproduce the non-commutative U{1) theory, we substitute the com- 
pleteness relation 



Ar2 

into the action (B.TI) and obtain 



-1 J]A(x) = I;v (4.12) 



Stek = -j;^J2Y. [ (f^r,f/trt + h.c.) A{x)] . (4.13) 

X fj, 

As in the context of twisted reduced models, the matrices act as lattice shift 

operator, and thus they behave as discrete derivatives e'*'^'". Indeed one can show 
from above that 

r^A{x)Tl = A{x-afi), (4.14) 
from which it follows that shifts on the fields are represented as 

U{x + afi) = ^tr(^) {T,UTlA{x)) . (4.15) 
Therefore, we can rewrite the action (B.TI) as 



STEK = -p^^[K{x)i.U*{x + afi)+U{x + afi)*U*{x) . (4.16) 

X fJ, 

4.1 The non-commutative theory 

The theory described in eq. (|4.16| ) is naturally U{1) left and right invariant, i.e., 
given a constant field g G U{1) the action is invariant under the transformations 
(where the ordinary product is intended) 

U{x)^g-U{x), U{x) -^U{x) ■ g. (4.17) 

Let us turn to the commutative continuum version of the theory, whose field we 
call u[x). First of all we perform the following substitution, dictated by the U{1) 
unitarity condition 

u{x)u*{x) = 1 u{x) = e''^(^\ ipeR. (4.18) 
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The action then reads, up to the couphng constant 



S 



d X d,. 



d^x d^(p{x)d'^(p{x), 



(4.19) 



and therefore the theory is equivalent to the theory of a free massless real field. 

Turning now to the non-commutative continuum theory of eq. ( p.l8| ), we notice 
that if the field lA (x) decreased sufficiently rapidly at infinity we could integrate by 
parts and turn the star product into the standard one. The action would thus be the 
same as in the commutative version. Naturally, the field is subject to the constraint 
of star U{1) unitarity, which reads 



U{x)i.U*{x) = 1 = 

(4.20) 

Uix)U*ix) + Zn=i ilT ^^'^'^ ■ ■ ■ ^^"'"^n ■ ■ ■ ^^Mx)^,, ■ ■ ■ d,JA*{x). 

This condition naturally implies that the theory is not a free theory as in the com- 
mutative case, although the action would be formally the same. Moreover it is not 
clear if a scalar complex field subject to the condition ( |4.20| ) can at the same time 
satisfy the rapidly decreasing condition needed to integrate the non- commutative 
action by parts and neglecting the boundary behavior. Therefore, we will still use as 
the reference action of the continuum theory we wish to study eq. ( p.l8|) . 



4.2 Physical interpretation: the brane vacuum 

The low energy effective action for a p-brane in the presence of nonzero constant B^y 
field along the brane is given by the dimensional reduction of the 10-dimensional 



non-commutative Yang-Mills model to the brane world volume [^. In particular, 
in the limit of large field B^y, i.e. when 

a'\\B^y\\ » \\g,yl (4.21) 

the non-commutativity parameter 6^^ is given by 

Qf,u ^ {^B-^Y'^. (4.22) 

In the case of a brane-antibrane pair of a non-BPS non-stable brane, one finds in the 
spectrum of the effective theory also tachyonic modes, described by non-commutative 
scalar fields T with tachyonic potential V{T) ||TT |. 



In the trivial non- commutative gauge field background, the part of the action of the 
brane system containing the tachyonic mode is given by the same action describing 
a non-commutative Higgs-like model of charged scalar fields, i.e.: 

S = I d^+^x (^9^0* 9^0^ - \/(0^0t) j , (4.23) 
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where V{-) is a potential with a nontrivial v.e.v.: |0p= some constant, and the field 



transform in the bi-fundamental representation of the U{1) gauge group . 

A point in the (true) vacuum of the field (j) can be parametrized by an element of 

the non-commutative U{1), 

(h ^ U -k (h , ^ 

where U and U'^ satisfy [/ ^ f/^^ = f/^^ * f/ = 1. 

To get the action describing the dynamics along the valley of the potential V in 
terms of the field [/, one has to take (f) constant at the minimum of the potential, 



perform the transformation ([4.24|) and declare U dynamical 0. The action for the 
Goldstone U- field is then 

S = ^j d^^Y'd^U * dM~\ (4.25) 

where 1 / = (p-k (p^ . 

Our model then corresponds to the case p = 1, with Euclidean metrics, and with 
the obvious identifications for the couplings A, P and for the constant field and 
the non-commutativity parameter of eq. ( [4.11| ) according to (|4.22| ). 



4.3 Double scaling limit 



From ([4.11|) we see that in order to take the continuum limit of the model in such 



a way that the dimensionful non-commutativity parameter 6 (which in the present 
case is just a real number) is fixed, we must fix the quantity a'^N. It is clear from 
above that we have to send N to infinity if we want a to go to zero and the dimension 
of the lattice to go to infinity. In order to set a — > we have to tune somehow the 
coupling constant (3. From renormalization group analysis of the beta function of 
chiral models it is known that 

a ~ A-^e""^, c = Svr. (4.26) 

It is questionable whether a relation like ( |4.26| ) is valid in the context of the TEK 
reduced model, because whether the equivalence between principal chiral models and 
TEK reduced models holds is in itself a nontrivial question Moreover equation 
( [4. 261 ) is strictly valid only in the planar limit. Nevertheless, we would like to propose 



to assume such a relation, and to numerically verify its consistency. 

An analogous assumption is made in [|l3l, where the relation between a and (3 is 



taken to be the known Gross- Witten planar result. The numerical results presented 
in 



13|, incidentally, strongly confirm the validity of such hypothesis. 



Whether a different value for c from the one indicated in ( [4.26| ) or a different 



functional dependence would lead to similar results to what we present in sec. ^ is a 
non trivial legitimate question. 
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What we propose to do is therefore to send and (3 to infinity in such a way 
that = N ■ e~^^'"^ is kept fixed. Numerical analysis indicates that indeed finite 
and (5 effects tend to compensate in such a limit, in a manner similar to the one 



obtained in |13] for the 2D EK model. 



4.4 Correlation Functions 

Typical objects that can be studied in numerical simulations are correlations func- 
tions. In particular, the easiest one to compute turns naturally out to be the two 
points function. In particular, given a lattice site n = (ni,n2), we define, in the 
reduced model, the following function, which is nothing but the translated version of 

G(n) = ^(Metr(^) (f/(n)t/(0)t)) (4.27) 



via the substitution of eq. ( p.3|) , namely 

GtekH ^ ^(Metr(^) (r^r^ ^(r|)"^(rl)"^f/t) ). (4.28) 



Again substituting the completeness relation ( [4.12[ ) we get 



GtekH = ^ J](Re(W(a; + n)^W*(x))). (4.29) 

X 

Since the non-commutative theory is defined on a lattice with A^^ sites, this expression 
defines the average of the two point function over all the possible lattice sites, and 
thus gives a coherent expression for the two point function of the theory. 

Incidentally, the internal energy of the model is given, up to constants, by 
Gtek{1,0) + Gtek{0A). 



5. Numerical Results 



What should we expect from numerical Monte Carlo analysis? First of all, if the 
double scaling limit we take is correct, we should expect a coherent superposition of 
the behavior of the correlation functions a.s N,i3 —>■ oo. If this is the case, it would 
strongly indicate the correctness of the procedure. Secondly what we should not 
expect is the typical behavior of the finite A^ corrections found in the reduced model 
in the strong coupling regime 0. In the double scaling limit one takes into account 
also the non-planar diagrams, and it is therefore not clear how the finite A^ and f3 
effects could manifest. 

We used a Metropolis algorithm to update the N x N matrix U. Trial matri- 
ces were selected by multiplying the actual matrix f/ by a random SU (2) matrix 
embedded in SU{N), choosing randomly among the A^(A^ — l)/2 SU{2) subgroups, 
and a U{1) random phase. More precisely, once the SU{2) subgroup and the U{1) 
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phase were randomly chosen, we performed ten Metropohs hits with an approximate 
acceptance of 50%. Each SU(2) updating requires 0{N) operations. The number 
of S't/(2)-subgroup updatings per run was 0(10''). We should also mention that for 
the values of N we investigated (see Tab. [l|), we observed some problem of thermal- 
ization when using completely random configurations (for example, constructed by 
multiplying A^(A^ — l)/2 completely random SU(2) matrix embedded in a x 
and associated to different subgroups, times a phase) as starting point of our simu- 
lations. In particular, we noted a worsening with increasing j3. So, as starting point 
we used either moderately random matrices or the unity matrix. It is well known 
that a simple Metropolis algorithm does not provide a particularly efficient method 
to simulate a statistical system. Our choice was essentially due to the fact that the 
reduced action is quadratic in the matrix variable U. Moreover, it does not lend 
itself to a linearization by introducing new matrix variables, as in the case of the 
reduced TEK gauge theory [1^ . 

We numerically studied the theory for 
two particular values of ^9. Our choice was 
motivated by two facts. First of all it seems 
unnecessary to test the existence of a dou- 
ble scaling limit in the strong coupling re- 
gion (i.e. for /3 < /3c ~ 0.3058) since the 
reduced theory is actually under control in 
that regime, and it has been shown ^ to 
accurately reproduce the standard, commu- 
tative parent theory. Secondly, we must ad- 
dress to sufficiently high values of A^, be- 
cause on the one hand the limiting proce- 
dure is expected to be sensible in the large 
limit, and on the other hand the dual lat- 
tice of the reduced theory has dimensions proportional to A^. Therefore, in order to 
avoid what we can legitimately call finite size effects, we had to resort to high values 
of A^ (the correlation length of the non reduced model is of the order of some lattice 
spacing in the region we investigated). 

We choose to take the values reported in Tab. |l], and to take A^ > 50. A reasonable 
statistics for our Metropolis algorithm limited the highest A^ value to A^ < 120. 





A^ 


/5 


^1 


= 9.93 


10-6 


60 


0.31063 








70 


0.31370 








80 


0.31635 








100 


0.32079 








120 


0.32442 




= 1.14 


10-6 


50 


0.35000 








80 


0.35935 








100 


0.36379 








120 


0.36742 



Table 1: Set of values of {N,[5) at fixed 
{} which we investigated numerically 



The correlation functions we studied was defined on the lattice as 
Gin) = . 



(5.1) 



In Fig. m and |] we show our numerical results. 

For large n, we see that asymptotically the obtained values for Q{n) agree within 
error-bars. We should mention that we numerically studied also the so called diagonal 
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Figure 1: Numerical results for Q{n) 

two points correlation function on the lattice, defined by 

^d(v^-n) = GTEK{n,n). 



(5.2) 



The results we obtained for Qdin) show an analogous asymptotic superposition within 
error-bars for large n in the double scaling limit, with values coherent with what was 
found for Q{n). In fact, the lattice definition of the two point correlation function 
should not depend on either of the two definition one takes, and this can be viewed 
as another confirmation of the validity of the outlined scheme. 

In all cases, what we find supports the expected numerical scenario, and con- 
firms that the procedure described above indeed yields a sensible nonperturbative 
definition of the non-commutative theory we described. 
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Figure 2: •d2'- Numerical results for Q{n) 
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